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ly^ • Abstract 

' In this paper, forward-backward asymmetry and differential cross section of top quark in 

flavor violating Z' model up to (^(a^ax) at Tevatron are calculated. In order to account 
^ _ for the top observed large forward-backward asymmetry, the new coupling gx among Z' 

and quarks will be not much less than strong coupling constant gs- After including the 
Oh' new higher order correction, the differential cross section can fit the data better than those 

, only including the leading contributions from Z' , while the forward-backward asymmetry 

is still in agreement with the measurement. 
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O : I. INTRODUCTION 

in 

I As the heaviest fermion in the standard model (SM), top quark is thought to be closely 

^ I related to the mechanism of electroweak symmetry breaking and physics beyond the SM 
(BSM). In the last two years, DO and CDF Collaboration have measured the forward- 
backward (F-B) asymmetry (Afr) of top quark at the Tevatron [H-I^. SM predictions 
have been estimated in Refs. |#^. In the SM, the asymmetry arises from the interference 
among virtual box and the leading diagrams for the process qq — t- tt, as well as the 



O 
O 



X 



^ ■ contributions from qq — )■ ttg. The present experimental measurements and SM theoretical 



predictions are listed in Table [T]in the lab {pp) frame and the center-of-mass (cm.) frame 
of the top quark pair (tt), respectively. From the table we can see that the CDF measured 
A^p^ is consistent with Ap^, if the theoretically expected dilution of 30% is included 
However the SM predictions is significantly smaller than the observations. 

TABLE I: A^ collection of experimental and theoretical results of Aps of top quark at the 
Tevatron \u- 
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DO (0.9fb-^) 


CDF (1.9fb-i) 


CDF (3.2fb~i) SM Theory 


Att 
^FB 


0.19 ± 0.09 ± 0.02 (exclusive 4 jet) 
0.12 ± 0.08 ± 0.01 (inclusive 4 jet) 


0.24 ±0.14 


0.078(9) 


APP 
^FB 




0.17 ±0.08 


0.193 ±0.065 ±0.024 0.051(6) 



Recently some theoretical progress has been made both in the SM and the BSM, in 
order to explain this novel signature. In the SM, soft gluon resummation effects [Tj] have 
been scrutinized. However, the prediction involving resummation effects does not change 
the asymmetry at ^(al) greatly Many BSM models, for instance, supersymmetry, 
extra dimension and left-right model have also been considered jln]. New particles such 



as exotic gluon G', extra W or Z' bosons and extra scalar S are introduced. All these 
new models should produce the required asymmetry while keep other observable qualities 
to be consistent with measurements. Among which, the ti invariant mass distribution is 
an important measurement to constrain the new models. In order to distinguish different 
models, as depicted in the Ref. [l2|, higher order effects in these new models are important. 

In this paper, we are interested in a BSM model named the flavor violating Z' model 
(FVZM) [lO(]. Observed asymmetry can be generated by introducing a right-handed 
coupling among the Z', the top and up quarks C 9 gxZ^'u^^^^-Y-t + H.c. A detailed 
analysis based on leading order (LO) contributions has been given in Ref. |lO|. For 
the suitable parameters, while the asymmetry can be generated, the tt invariant mass 
distribution does not fit the observation well. Therefore it is quite interesting to analyze 
the asymmetry and tt invariant mass distribution after including higher order effects. 

The paper is organized as follows. In Sec. HIl Born, virtual and real corrections are 
calculated analytically till to ^{a'^ax)- In Sec. IIIIl numerical results for differential cross 
sections and forward-backward asymmetry as a function of the top quark pair invariant 
mass Mft are presented and compared to the experimental data. In Sec. llVt we give a 
short conclusions and discussions. 



II. ANALYTICAL CALCULATION UP TO ^{a'^ax) 

In this section, we will present the analytical formula to calculate the top forward- 
backward asymmetry, as well as the differential cross section up to ff{a1ax)- The cor- 
responding Feynman diagrams of subprocesses up to i^^a^ax) in FVZM are depicted in 
Figs. 1-5. In order to account for the top large asymmetry, the new coupling gx will be 
not much less than Qg. Thus the relevant amplitude for ti final states can be written, in 
perturbation series of couplings, as 

= fsOs + fxax + fla] + flx^sax + ■■■ (1) 

with as = (7^/(47r), ax = 5'l:/(4vr), and /'s the corresponding form factors. Squaring the 
amplitude we obtain 

|^*f = \fs\'al + 2Mif:fx)asax + \fx\'al 

{f:fl) al + 2^ {rjl^ + r^fl) alax + ■■■ (2) 
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In order to cancel the infrared divergences, the corresponding gluon radiation processes 
should be included. The amplitude for tig finals states can be written similarly as 



(3) 



Squaring this amplitude we obtain 



(4) 



In the SM, the asymmetry arises from the ^{ds) term. In the FVZM, new contributions 



till to ff{ax) are calculated in 10|]. In this paper, the extra contributions at ff{alax) 
will be calculated. The SM uu — )■ ti, dd — )■ ti up to QCD NLO and gg — )• ti up to QCD 
LO contributions are recalculated though their analytical expressions are not shown in 
this paper. 



A. Contributions up to & {ctx) 

Typical Feynman diagrams, which contribute to the amphtude up to 0" (a^), are shown 
in Fig. [H 



t 




u 



FIG. 1: Typical Feynman diagrams with contributions for form factors fs and fx- 



The form factors of 2^ [f* fx) and with spin- and color-summed (same for the 
following form factors) are given by 



(/:/- 



X) 



\fx\' 



1^""^^% , \m\ + (2m|, + s-2t) m\ 
+ {t^ - 2m|,(s + 2t)) m\ + 2m|,(s + tf \ , 



m. 



m%, {t — 
-8m|,(s + t)mf + 4m|,(s + 



(5) 



(6) 



where Ca = 3, Cp = 4/3 and s = {pi+ P2Y, t = {pi — ki)^ are the Mandelstam variables. 



B. Contributions at ^ (a'^ax) 

The corresponding Feynman diagrams related to {a'^ax) are shown in Figs. 2-5. In 
order to regulate the divergences, dimensional regularization is adopted with D = A — 2e. 
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Infrared (IR) and ultra violet (UV) divergences are represented by 1/ ejR and 1 / euv respec- 
tively. The wave function renormalization constants are determined by the on-mass-shell 
scheme while MS scheme is chosen for the strong coupling constants renormalization. 
The calculations are carried out with the help of FeynCalc 18|], FormCalc 19| and QCD- 
loop [2^]. At hadron collider, in order to eliminate the collinear singularity, factorization 
should be carried out. In this paper MS factorization is adopted, as shown explicitly 
below. 
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(b) 

FIG. 2: Typical Feynman diagrams with contributions to form factors f^-^. 



Once the renormalization is carried out by adding appropriate counterterms, the am- 
plitude is free of UV divergences. The counter terms are given explicitly in Appendix 
A. However there are still IR divergences in f^x and f^. In order to eliminate the re- 
maining infrared divergence, the gluon emission processes should be included as shown in 
Figs. 4 and 5. Two cutoff phase space slicing method {21J is applied for such processes. 
First, small parameter Ss is introduced to separate the final phase space into soft part and 
hard part. Second, small parameter 6c is introduced to divide the hard part into a hard 
collinear part and a hard non-collinear part. These three parts are calculated separately. 
Their summation should be independent of the two small parameters Sg and Sc- 

The remaining IR divergences appear in 2^(/*/j^ + fxfl)(^1(^x [cf. Eq. [2]] can be 
expressed as 

'^^Ulflx + fxflKax = Vf + V,— + V2^, (7) 

e/fl (^IR 












u 8 t 



u S t 



u g 



FIG. 3: Typical Feynman diagrams with contributions to form factors fg. 
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FIG. 4: Typical Feynman diagrams with contributions to form factors /, 
where Vf indicates the finite part and IR coefiicients Vi, V2 are 

+91og(^) + lOlog(^) - 161og((^) -2logitlI!^)} 

37r2 r(l-e) s{t-ml,) ' 



V2 = [2^(/;/x)'a.ax] 



-Was (47r)" 

SttCaCf r(i-e) • 







FIG. 5: Typical Feynman diagrams with contributions to form factors fx- 
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where f3 = yl — Ami/ s ^i^nd /x is an energy scale introduced in dimensional regularization. 
Here 

is different from 2^{f*fx) in Eq. [5] as Goldstone boson contribution is ignored. 

The IR poles in Eq. [7] have two physical origins: soft & collinear divergences. The 
double pole indicates an overlap between soft and collinear divergences, and the 

divergences can be eliminated by including contributions from the soft region for gluon 
emission processes and parton distribution function (PDF) redefinition 2l|]. Gluon emis- 
sion process in soft region can be calculated by the eikonal approximation method (the 
details are given in Appendix B), which can be expressed as 

2M UTfx) =Sf + R,— + (10) 



where Sf indicates the finite part and IR coefficients R2 are 

log(i±f ) + 16 + 181og(^) - 32 log 5, 
+201og(^) - 321og(M^) - 41og((^±^)} 

where 2^{f*fx)^ is a D-dimension version of 2^(/*/x)', 



^D_.oo_2^ ^ (1 - e/i;)K + + - ^ti^ + 2t) + sjt - m1)ern] 

s{t — m|/) 



PDF redefinition in soft region can be written as [2l| , 

[2^(/:/x)V«x]^:p^^(3C^ + 4C^log5,)(log(4) - — )■ (13) 
47rr(l-e) /i^ ejR 

IR cancelation is realized by adding Eqs. [7] and [10], and subtract 2 times of expression 
in Eq. [I3]. 

For the hard collinear part for gluon emission processes, there are collinear divergences 
which can be eliminated by PDF redefinition in this region. The remaining finite part is 
in the form of a convolution integral with a splitting function. The detail of this procedure 



is described in Ref. 21 



The remaining hard noncoUinear part of gluon emission processes is finite and the 
integration is performed with a standard three body Monte Carlo code. 

During the calculation, there are three scales in the hadron level cross section, /i is 
an energy scale introduced in dimensional regularization. fij. is a renormalization scale 
introduced in the MS renormalization of the coupling constant, /i/ is a factorization 
scale introduced in MS factorization. The fi always comes with the divergences and is 
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canceled completely when the corresponding poles are subtracted. The dependence of our 
final results on /z^ & /i/ will be shown in the numerical results. 

The independence of the total cross section with 6s and 6c has been checked in the 



situation 6c << 6s, as suggested by the two cut off phase space slicing method in Ref. [21 



III. NUMERICAL RESULTS 



In this section, we will present the numerical results and compare them with the 
experimental measurements. We choose cteq61 for leading order calculation and cteq6m 
for higher order calculations. The scales /x^ and Hf are set to be equal and ^^(m^) = 0.118. 
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FIG. 6: Differential cross sections as a function of M^j with fi^ = fJ-f = fnt- Here "QCD Born" 
and "QCD NLO" represent the results in the SM at leading order and next-to-leading order in 
QCD. "QCD Born + Z' Born" and "QCD NLO + Z' NLO" stand for the predictions in FVZM 
up to {ox) and {a^ax) respectively [c.f. Eqs. [2]and|l]]. 



Differential cross sections as a function of Mu are shown in Fig. |6l Histograms are 
drawn here in order to compare conveniently with the experimental measurements jsf. 
The parameters in the FVZM are taken to be ax = 0.024, Mz' = IGOGeV which is the 
best point [10] to account for the top asymmetry. From the figure it is obvious that 
the NLO QCD prediction is in good agreement with the data except the bin around 400 
GeV. It should be noted that even the multiple soft gluon radiation effects are included, 
the discrepancy remains (c.f. Ref. [3]). However the top quark asymmetry at NLO 
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QCD is much less than the measurement. After including the contributions from the 
leading diagrams from Z', for the favorable parameters, the top quark asymmetry can 
be generated. However the differential cross section does not agree with measurement 
well. From the figure, we can see the prediction is lower than measurement for small Mti 
region while higher for large Mti- After including the higher order effects, the differential 
distribution will be better while the top asymmetry can be generated. These behavior can 
be understood as follows. In the vicinity of ti production threshold region, the significant 
contributions comes from the interference among QCD and extra Z' Feynman diagrams. 
Such contributions will decrease the cross section. In the higher Mjf region, the square of 
Z' diagrams become significant and they will uplift the cross section. The deviation from 
NLO QCD prediction will be soften after including the higher order contributions. 




FIG. 7: Differential cross sections da/dMii as a function of M^i with /x,. = fif = 0.5'mt,mt,2mt, 
respectively, for 4 sets of parameters. Other conventions are the same with Fig. [6l The factor 
k = a^^'-'/a^^ is also indicated. 

In Fig.[7]we show the differential cross sections da/dMu with = l^f = 0-5mt, rrit, 2mt 
respectively, for four sets of typical parameters, namely different Mz' and ax- From the 
histograms, we can see that the above-mentioned improvement after including the higher 
order effects is universal. We calculated the k factor, which is defined as = a^^^ /a^^ . 
Here a'"^ and a^^^ are the cross sections up to 0'{a\) and ff{agax) respectively. For 
four sets of parameters, the k-factor is equal to 1.03, 0.996, 0.998, 0.992, respectively. 
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for jjLr = jJif = rrit. It is obvious that the NLO contributions mainly change the shape of 
distribution. As for the scale dependence, the results of "QCD NLO" and "QCD NLO+Z' 
NLO" are about the same size. "QCD NLO+ Z' Born" is significantly smaller than them. 
"QCD NLO + Z' Born" scale dependence is small because QCD NLO and Z' Born have 
opposite fir/fJ'f dependence. 




FIG. 8: Forward-backward asymmetry distributions as a function of M^j with fXr = = mt. 

Histograms of the forward backward asymmetry A^p^ as a function of Mti in ti rest 
frame are drawn in Fig. [HI where A^p^ in each bin is defined as 

- iv(Ay > 0) - iv(Ay < 0) 

N{/\Y > 0) + N{/\Y < 0) ^ ' 

where AY = Yt — Yi denotes the difference between the t and i rapidities. The total A^p^ 
is calculated to be 15.8%, by summing A^p^ in each bin multiplied by their corresponding 



weights. To compare directly with experimental data [22|, we also draw the so called 
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"above" and "below" ApB distribution at the bottom of Fig. [8] in which ApB measured 
or calculated for Mu above or below a certain value. It should be noted that the experi- 
mental data are measured in the pp lab frame. As the simplest approximation, we utilize 
the relation A^p^ ~ O.TA^p^- Obviously more measurements are needed to decrease the 
experimental uncertainties in order to confirm/exclude the FVZM. 

Total cross sections and total A^p^ and are shown together in Table [III For total cross 
section, the Z'-born contribution decreases the NLO QCD cross section. Including the 
Z'-NLO corrections makes the cross section even smaller although these corrections are 
not significant. On the contrary is sensitive to Z'-NLO correction and can drop 

about 30% from the Z'-Born value. 

TABLE II: A^p^ and total cross sections with mz' = IGOGeV, ax = 0.024, /i,. = Hf = m^. 



QCD NLO QCD NLO+Z' Born QCD NLO+Z' NLO 
6.8 22.2 15.8 

Total cross section(pb) 6.29 5.52 5.13 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper, we calculate the top quark differential cross section and asymmetry 
up to ^(a^ax) in a flavor violating Z' model (FVZM). In the FVZM, the leading Z' 
contribution can induce the measured top asymmetry, while the differential distribution 
of Mft does not fit measurement well. After including the higher order contribution, the 
differential distribution can be improved while the top asymmetry is still in agreement 
with the observed value. 

QCD soft gluon resummation effects for the top quark pair production in the SM have 
been considered in Ref. Such effects do not change the M^f distribution significantly. It 
is expected that resummation effect in the FVZM is similar to that in the SM because the 
internal Z' contributions has nothing to do with the soft gluon radiations from external 



quark legs. Such effects are under investigation [23 
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Appendix A. Renormalization Constants 

Renormalization constants are needed when calculating Fig.|2]and Fig. [31 6Zu, SZt, 6Za 
corresponds to up quark, top quark and gluon on-mass-shell renormalization constants 
respectively. 6Zg is the coupling renormalization constants. 
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^Za = trtS)((i^^ - i^^^'/ - l^FUHf)^ - {ICa - ITfu.j)^^ (15) 
+|THlog(f) + log(^)+log(^))), 

where = 4/3, Ca = 3, Tp = 1/2, ri// = 3 is the number of hght quark flavors, Uhf = 3 
is the number of heavy quark flavors, n is the energy scale introduced in dimensional 
regularization, /z^ is the renormalization scale. 

There are new contributions to top and up quark fleld renormalization constants, when 
we calculate Fig. |2l Top quark fleld renormalization constants SZY , SZf are calculated 
from Z' induced top Self-energy — (13], as showed in Fig. [9l 



z t 




II t 



FIG. 9: Self energy diagrams for counterterm calculation. 

5ZY = Sr(p^)|p2=„2 + 2m?^(Sr(p2)|,2=„. + Sf (16) 

5Z^= S^(p2)|^,^^,. 
The counterterm for tig vertex is written as 

(-2(7.TV)(5^r + ^^*V), (17) 
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where the vector and axial vector parts are 



xyv _ _JM)^_gxL- 

"^t - r(l-e) 32m^7r2 



"^t — r(l-e) 32mt27r2 



+ [-mt^ + (m,.2 - ml)Bo{0, ml m,,^) 
+ {mt^ + ml- m,>^)Bo{mt^, ml, m^^)]}. 

(18) 

Field renormalization constants of the up quark can be calculated similarly and the 
counterterm for uug vertex is written as 

{-zgsT^Ym: + SZ^j'), (19) 

where 
^yv — __(Ml_22f!_ 

"^M r(l-e) 32m27r2 



6Z. 



+ {mt'^ - m^/2)5o(0,mt2,m2/2) + (-mj^ + ml + m,,^)Bo{ml,mt'^ ,m;,,'^)]}, 

(20) 

_ (4^)^ ax^ 

u T{l~e) 32m27r2 

+ [-ml + (m,,2 _ mt^)Bo{0, mtl m,,^) 
+ {mt^ + ml - m,,^)Bo{ml,mt^,m,''^)]}. 

Appendix B. Soft part of the real gluon emission cross section 

Soft real squared amplitude |A^q^{^jf^P, which is expressed as [2^ {fl* fx) '^s'^^lsoft 
Eq. [To], can be obtained by the interference of diagrams in Figs. [Hand [5] with requirement 
that the gluon's energy is smaller than 5s\fs/2 . |-A^!^{^jfoP can be expressed as 



qq-^ttg I 

4 



= l-M,.--d\E (21) 



where |A^gg_>.tip = 2^{f*fx) OigOix is the interference term of the two born diagrams 
in Fig. [U Cij = Cji is the color factor of interference terms with one gluon emitting from 
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external leg i of one diagram and from external leg j of another diagram. 

Cl2 = C*14 = C*23 = C*34 = ~Cp/2, 



(22) 

11 = ^22 = '-'IS = "--24 = "--33 = '-'44 = '-^A'-^F- 



/^2 

L/ii — L/22 — (--IS — L'24 — '-'33 — '-'44 — ^ 



Co = CaCf is the color factor of the interference of the two diagrams in Fig. [H Sij = Sji 
are the soft factors of the corresponding interference terms. They are calculated by using 



eikonal approximation method. According to Ref. [21 



C _ as (47r)" /fj.^ 
"^i-i 2-K r(l-e) V s 



■K^ S 



(23) 

where r]i is a sign which is positive for outgoing fermion or incoming antifermion, and is 
negative for incoming fermion or outgoing antifermion. 

5*11 = 5*22 = 0, 

^12= tT^(4-21og(^)^ + 21og^(^)-^), 

^13 = S2. = trtSyllt - + 

-ilog^(fil) + log^(^) + ilog^(^) 

+ log(^) log( (i^^) + U2i^^^) - l^2i-^^^) - S), 

= S2S = t rtSy{-|4 + [log(^) + I log(^^^)]^ (24) 

- log(^) log( (^^f^) - h2{^^-^) + U-^^^) + fj}, 
^33 = ^44= trtSy(^ + ^log(f^)-21og(-f^)), 

^34=-^{-^l0g(|i^)^ 

-^[ilog^(fS)-log(^)log(|^) + /^2(|fr)]}^ 
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in which (3 = yl — ^mf / s and 9 is the angle between the incoming u and outgoing t 
quark in tt rest frame. 
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